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Abstract. In this paper, we discuss a two-dimensional numerical study of isothermal high Mach number Bondi- 
Hoyle-Lyttleton flow around a small accretor. The flow is found to be unstable at high Mach numbers, with the 
instability appearing even for a larger accretor. The instability appears to be the unstable radial mode of the 
accretion column predicted by earlier analytic work. 
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1. Introduction 

The Bondi-Hoyle-Lyttleton (BHL) accretion problem 
concerns the accretion rate of a point mass moving at 
constant velocity through a u niform gas cloud. Th is prob- 
lem was first considered by iHovle fc LvttletonI (^939). 
Considering Newto nian trajectories in a pressureless fluid, 
iHovle fc LvttletonI concluded that material with an im- 
pact parameter {Q less than 



C<C: 



2GM 



HL 



(1) 



would accrete onto the point mass, M {v^ is the velocity 
of the gas relative to the point mass at infinity). If the 
density of the gas is (at infinity) then the accretion 
rate onto the point mass will be 



A/, 



HL 



(2) 



For a full discussion of the probl em, see lEdgail lj2004 
Numerical studies have shown fsee lEdgail l|2004|) . and ref- 
erences therein) that, although the accretion rate pre- 
dicted by equation [3 is fairly accurate, the flow geome- 
try is somewhat more complicated, due to the effects of 
pressure. 

Bondi-Hoyle-Lyttleton accretion is of particular in- 
terest in star formation studies, since protostars in clus- 
ters are thought to undergo competitive accretion, at 
rates similar to that given by equation [3 ((Bonnell et al.l 
However, under such conditions, is low, and 
so CjjL is very large compared to the radius of the ac- 
cretor (i?acc)- This region of parameter space is poorly 
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studied, due to the larg e computational expense involved. 

observed a flow instability while 
studing the effect of radiative feedback on Bondi-Hoyle- 
Lyttleton flow, and this paper looks more closely at this 
instability. The computational method is discussed in sec- 
tion|3and the results in section^ Finally, the conclusions 
are presented in section^ 

2. Method 



The Zeus-2d code of IStone fc NormanI l)l992(l was used 
to perform the calculations presented here. The author 
parallelised the code using OpenMP, and found that it 
scaled well to up to eight processors. Extra routines to 
monitor accretion rates and drag were added to the main 
code, and calculations were performed in spherical polar 
(r, 0) co-ordinates. 

A linear grid of 100 cells was used in the 9 direction, 
while the radial grid was logarithmically spaced. The size 
of the innermost radial cells was chosen to make the in- 
ner portions of the grid as 'square' as possible (to avoid 
possible truncation errors). For the inner boundary con- 
dition, three alternatives were tried: Zeus' own 'outflow' 
boundary, a constant (small) density, and the density set 
to a small fraction of the innermost 'active' cell. The outer 
boundary lay at IOChl (^^ defined by equation^ with a 
uniform inflow imposed on the upstream side, and outflow 
conditions on the downstream side. 

The grid was initialised to a uniform flow, and the ac- 
cretor introduced over a number of timesteps at the start 
of the run. The simulations were run for as long as possi- 
ble, to ensure that any 'switch-on' effects died away. The 
goal was a crossing time of the whole grid, but this was 
not possible for the runs with the small accretor. 



2 



R. G. Edgar: High A4 BHL Flow around a small accretor 



1 .A r 




10 2G 30 40 

Simjiation tin^e 



Fig. 1. Accretion rates as a function of time for low Mach 
number flow around a large accretor 

In the simulations, dimensionless units were used. The 
incoming material had = 1, and the sound speed was 
set according to the required Mach number {Ai). The 
point mass was set to M = 1 (as was the gravitational 
constant, G). Using equation ^ we see that Chl = 2 for 
these choices. A value of = 0.01 was chosen (arbitrar- 
ily, since the flow was isothermal and not self-gravitating), 
giving Mhl = 0.126 in these units. 

3. Results 

Three sets of simulations were performed. For each of the 
sets, one run for each of the inner boundary conditions was 
made. First, a low Mach number {A4 = 3) flow around a 
large (i^a^c = O-OO^^hl) accretor was used to test the code. 
The second set used the same size of accretor, but with 
the Mach number increased to = 20. For the final set, 
the accretor size was reduced to O.OOOS^jjl. 

We will now discuss the results of each of these sets in 
turn. 

3.1. Low Mach Number, Large Accretor 

The accretion rate as a function of time is presented in 
figure^ After an initial turn-on period, the accretion rate 
settles to a fairly steady value, slightly larger than that 
predicted by equation |21 Changing the inner boundary 
condition does not affect the results significantly, as is 
shown in figure 13 The different boundary conditions are 
giving slightly different curves, but are very close together. 

The gravitational drag force on the accretor for this 
set of runs is plotted in figure 13 This is plotted in imits of 
-^HL^oo which is the natural scaling for the problem. The 
drag force rises smoothly to around ISM^^v^. 

Inspection of the density field shows a wake with a 
fairly broad opening angle, but no bow shock - as is typical 
for isothermal BHL flow. None of these results ar e partic- 
ularly surprising, when compared to the work of iRAifferil 
l)l99(j) (model 'GS' is closest to that presented here). From 
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Fig. 2. Detail of figuren]to show differences due to bound- 
ary conditions 
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Fig. 3. Gravitational drag force in units of M-^j^v^ for 
low Mach number flow around a large accretor 



this set of runs, we can conclude that the code is behaving 
well. 

3.2. High Mach Number, Large Accretor 

The accretion rates for the three runs of high Mach num- 
ber flow around a large accretor are presented in figure 21 
The behaviour is very different to that shown in figure ^ 
with a large oscillating accretion rate. We can see that 
changing the boundary condition does change the flow his- 
tory, but the change does not appear to be significant - 
the instability is similar in all cases. 

Figure |S1 shows the corresponding gravitational drag 
for this set of runs. As might be expected from the accre- 
tion rate curves, it is quite noisy. Again, there are slight 
differences due to the inner boundary condition, but these 
do not seem to be significant compared to the general in- 
stability of the flow. 
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Fig. 5. Gravitational drag (in units of Mhl''^oo) a func- 
tion of time for high Mach number flow around a large 
accretor 



The wake's opening angle is much less than for the low 
Mach number case, as one would expect. The oscillating 
accretion rate appears to be the 'ov ersta bility' of rad ial os- 
cillations described bv lCowiel l|l977ll and lSokejl(l99(t) . The 
analysis was performed for a pressure-free (ballistic) flow, 
and the high Mach number isothermal nature of th ese cal- 
culations is a close appr oximation of this condi tion. 
extend ed the analysis of lBondi fc Hovld l)l944) (see 
l(2004l) for a slightly more detailed exposition) to time de- 
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pendent flow. Perturbing the steady state solution and ap- 
plying the WKB approximation (that is, the wavelength of 
the perturbation is assu med to b e short compared to the 
scale on which it varies) . FCowid fou nd that the accretion 
column should be unstable. Soker showed that viscosity 
could suppress this instability (which has not happened in 
the present work) . Howe ver, the y disagreed on the physical 
nature of the instab ility: ICowid attributed it to the incom- 
ing material. ISokerl asserted that the equation showing the 
possibilty for instability (the dispersion relation from the 
WKB approximation) did not contain terms relating to 
the inflowing material, and hence the instability was in- 



trinsic to the flow along the line. Settling this dispute is 
beyond the scope of the present work. 

Looking at a single accretion rate history curve from 
figure^ around two to three oscillations occur per simu- 
lation time unit. This is in accordance with the prediction 
of Cowie, who pointed out that the instability must grow 
faster than the characteristic accretion time (Chl/^oo = ^ 
in this paper), or it will be swallowed by the accretor be- 
fore it can grow. Examining the velocity structure of the 
wake (figure El for a sample curve) , oscillations in v^. are 
ob served superimposed on the freefall curve, as predicted 
bv lCowi^ . The oscillation wavelength does not vary much, 
as required by the WKB approximation. The stagnation 
point is close to r = 3.5, which is unsurprising, given that 
the accretion rates are slightly higher than Mjj^ (which 
would place the stagnation point at r = 2 - cf Edgail 

iool). 

Unfortunately, computation of the power spectrum of 
the accretion rates (figured) is rather unenlightening due 
to noise. The spectrum peaks at a frequency of around 2 
to 3 (it is too noisy to be more accurate), and has a power 
law fall off at higher frequencies. There may also be a side 
peak at a frequency of around 0.5, but this is only cov- 
ered by a few points, and hence may just be noise. The 
very intense bursts of accretion seen in figure 01 seem to 
be associated with larger 'blobs' of matter being accreted. 
However, these are relatively rare, with only a couple ap- 
pearing over the entire duration of the simulation (and 
one of those was probably a 'switch-on' transient). 

3.3. High Mach Number, Small Accretor 

Finally, we come to the simulations of high Mach num- 
ber flow around a small (0.0005Chl) accretor. These sim- 
ulations were painfully expensive computationally. Each 
simulation in the flrst two sets required only about six 
days of wall time (on two processors). These final three 
each occupied two processors for eight months, but still 
did not achieve the same elapsed simulation time. As the 
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Fig. 7. Sample accretion rate power spectrum for high 
Mach number flow around a large accretor 




Fig. 8. Accretion rates as a function of time for high Mach 
number flow around a small accretor 

grid is extended towards the point mass, the gravitational 
acceleration gets stronger. There are more cells, the veloc- 
ities are higher, and the cell size is smaller. These three 
effects combine to cause the consumption of prodigious 
computing resources. 

The accretion rate onto the point mass as a function 
of time is plotted in figure |H1 The behaviour is generally 
similar to figure 01 with large fluctuations once the flow 
becomes established. There are differences in detail be- 
tween the three boundary conditions, but the qualitative 
behaviour is obviously the same. Figure El shows the grav- 
itational drag for this set of runs. Note that the gravita- 
tional force reaches larger values - as one might expect 
given the smaller size of the point mass (cf the Coulomb 
logarithm which app ears in th e calculation of dynamical 
friction. See e.g. Bin nev fc Trem ainc (1987)). If the flow is 
dense enough to make Mjjl interestingly large, then this 
drag force will alter significantly. However, the general 
pattern is similar to figure [SI showing fluctuations compa- 
rable to those of the accretion rate. 

The mean accretion rates for the two sets of simula- 
tions with high Mach number flow are presented in ta- 




Fig. 9. Gravitational drag (in units of M^-^Vr^) as a func- 
tion of time for high Mach number flow around a small 
accretor 

Table 1. Mean Accretion Rates in units of A^hl fo'' high 
Mach number flow for different inner boundary conditions 



Accretor Size 


Zeus 


Const, p 


Const, p Frac. 


Large 


1.47 


1.50 


1.48 


Small 


1.46 


1.46 


1.47 



bleQ] Shrinking the accretor does not seem to have made 
any significant difference to the results (beyond increasing 
the required computing resources by well over an order of 
magnitude). 

4. Conclusions 

Based on the results presented above, we see that high 
Mach number Bondi-Hoyle-Lyttleton flow around a small 
accretor is unstable. Although the average accretion 
rate is generally s i milar to that originally predicted by 
iHovle fc LvttletonI l)l939l) . short increases to rates several 
times this are possible. The instability seems to be related 
to the high Mach number of the flow - at lower Mach 
numbers, pressure must have a stabilising effect. So far as 
these simulations can determine, the instability is the ra- 
dial os cillat ion of t he acc retion column described bv lCowiel 
lll977l) and lSokeJ il99(t . 

Observationally, the instability may be visible for 
deeply embedded young protostars. For a 1 Mq protostar 
orbiting in a protocluster, the time unit in the present 
paper would correspond to approximately one year. The 
fluctuations in the accretion rate would correspond to a 
brightening of a factor of a few over a period of days to 
weeks, and dimming over a similar timescale. This is quite 
different from the FU Orionis phe nomenon which i s at- 
tributed to accretion disc activity fcf lBell et alJ l)200Cl|) and 
references therein) , where the outbursts rise on a timescale 
of months and last for years. 

In the author's opinion, there is little to be gained 
by the further study of the pure Bondi-Hoyle-Lyttleton 
problem. Shrinking the accretor size revealed no new be- 
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haviour, at a tremendous cost in computer time. Instead, 
the Bondi-Hoyle-Lyttleton model should be retained as a 
useful 'benchmark' and work focused on improving simu- 
lations of more realistic scenarios. 
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